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Abstract 

We study systems of Kuramoto oscillators, driven by one pacemaker, on d-dimensional regular 
topologies like linear chains, rings, hypercubic lattices and Cayley-trees. For the special cases of 
next-neighbor and infinite-range interactions , we derive the analytical expressions for the common 
frequency in the case of phase-locked motion and for the critical frequency of the pacemaker, 
placed at an arbitrary position on the lattice, so that above the critical frequency no phase-locked 
motion is possible. These expressions depend on the number of oscillators, the type of coupling, 
the coupling strength, and the range of interactions. In particular we show that the mere change 
in topology from an open chain with free boundary conditions to a ring induces synchronization 
for a certain range of pacemaker frequencies and couplings, keeping the other parameters fixed. 
We also study numerically the phase evolution above the critical eigenfrequency of the pacemaker 
for arbitrary interaction ranges and find some interesting remnants to phase-locked motion below 
the critical frequency. 

PACS numbers: 05.45.Xt, 05.70.Fh 
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I. INTRODUCTION 



Synchronization is an ubiquitous phenomenon, found in a variety of natural systems 
hke fireflies, where the males flash light in synchrony to attract the females chirping 
crickets P], or like neural systems, but also utilized for artificial systems of information 
science in order to enable a well-coordinated behavior in time. As an important special 
case, the coordinated behavior refers to similar or even identical units like oscillators that 
are individually characterized by their phases and amplitudes. A further reduction in the 
description was proposed by Kuramoto for an ensemble of oscillators js^, after Winfree 
had started with a first model of coupled oscillators Within a perturbative approach 
Kuramoto showed that for any system of weakly coupled and nearly identical limit-cycle 
oscillators, the long-term dynamics is described by differential equations just for the phases 
(fi (not for the amplitudes) with mutual interactions, depending on the phase differences 
in a bounded form. Kuramoto solved this model for the special case of all-to-all coupled 
oscillators, for identical interaction functions, and for eigenfrequencies, taken out of a narrow 
distribution. It is, however, this model, later named after him, that nowadays plays the role 
of a paradigm for weakly and continuously interacting oscillators. In a large number of 
succeeding publications the original Kuramoto model was generalized in various directions. 



for a recent review see jsl . In particular the eigenfrequencies were specialized in a way that 
one oscillator plays the role of a pacemaker with eigenfrequency different from zero, while all 
others have eigenfrequency zero 0]. In j^] the pacemaker was placed at a specific position, 
the beginning of a linear chain. Moreover the interaction range was changed from all-to-all 
to next neighbor interactions 0] ■ In [3] the system was studied numerically, and an abrupt 
change from long-to short-term interaction was suggested. 

In this paper we consider a system of Kuramoto oscillators, coupled on various regular 
lattice topologies, and driven by a pacemaker, placed at an arbitrary site of the lattice. 
The interaction range is neither restricted to next-neighbor nor to all-to-all couplings, but 
varied by a parameter a that tunes the decay with the distance between the interacting 
units. We analytically derive the common frequency Q of phase-locked motion and the 
upper bound on the absolute value of the ratio of the pacemaker's frequency to the coupling 
strength \uJs/K\ for next-neighbor and all-to-all couplings in case of (i-dimensional regular 
lattices with open or closed boundaries, in particular for a ring and an open chain. The 
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results show the dependence on the size of the system (total number of oscillators), the 
position s of the pacemaker (for open boundary conditions) and the (a) symmetric treatment 
of the pacemaker's interaction with the other Kuramoto oscillators (tuned by the parameter 
e^). The intermediate interaction range (0 < a < oo) as well as the phase evolution in the 
desynchronized phase (above the threshold for phase-locked motion) are treated numerically. 

The paper is organized as follows. In section II we define the model and identify the 
various special cases. Section IIIII is devoted to the study of phase-locked motion. We 
calculate the common frequency and the critical threshold for phase-locked motion for a 
linear chain and a —>■ oo fsection flll A|) . for a ring and a oo in section ITlI Bl the case of all- 
to-all couplings {a = 0) in section IIII CI Our main result about a switch to synchronization 
via the closure of a chain is stated in section IIIIDI Section IIII El presents our results of 
numerical simulations for the intermediate interaction range (0 < a < oo). In section IlII Fl 
we generalize the previous results to regular lattice topologies in higher dimensions, i.e. to 
Cayley-trees and d-dimensional hypercubic lattices). In section HVl we report on numerical 
results for the phase portrait above the critical threshold, in the desynchronized phase, in 
which the phases appear as a superposition of a fast motion, determined by the pacemaker, 
and a slow motion, determined by the interaction of all other oscillators. Finally we give the 
summary and conclusions in section |3 Details of the analytical calculations are postponed 
to the Appendix. 



II. THE MODEL 



The system is defined on a ^-dimensional lattice, where each site is a limit-cycle oscillator. 
The phase of the i-th oscillator obeys the dynamics determined by 

0. = kscos + (1 + S^.ses) - V ,yt = 0,...,N , (1) 

where 

Vi--=Y1''7^ >Vz = 0,...,Ar . (2) 

In particular (pi is the phase of the i-th oscillator, the index s labels the pacemaker with 
eigenfrequency Us 7^ 0, while all other oscillators have eigenfrequency zero. The parameter 
— 1 < < serves to tune the interaction of the pacemaker with the other oscillators. For 
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eg = the pacemaker is on the same footing as the other oscillators, being only distinct by its 
eigenfrequency uJs- For = — 1 its interaction is asymmetric in the sense that the pacemaker 
influences the other oscillators, but not vice versa (the pacemaker acts as an external force), 
so that it is on a different level in the hierarchy of the system. (In natural systems both 
extreme cases as well as intermediate couplings may be realized: a pacemaker being one out 
of an ensemble of oscillators or being distinct by additional features to the eigenfrequency.) 
The constant K > parameterizes the coupling strength, rji represents a space-dependent 
normalization, it replaces the factor that is used for all-to-all interactions and turns out 
to be essential for the a-dependence of the synchronization threshold, cf. section ITlI El The 
phases of the i-th and j-th oscillators interact via the sine function (as originally proposed 
by Kuramoto), of which two properties enter the analytical predictions: the sine is bounded 
and an odd function. The distance r^j, given in units of the edges of the lattice, enters to 
the power —a, a = gives the space- independent all-to-all coupling, originally considered 
by Kuramoto, a oo is understood to restrict the interactions to nearest neighbors. The 
label of the oscillator simultaneously specifies the location in space, see e.g. Figure [T] for a 
one-dimensional chain, where coordinate s specifies the position of the pacemaker. 

12 s-1 s s+1 N-1 N 

CHDO CH>0----0-0 

Figure 1: One-dimensional lattice with (A'^ + 1) Kuramoto oscillators, labelled by their coordinates. 

System (PJ contains the original model of Kuramoto as special case for = 0, a = 
and absent 5j (all oscillators with an eigenfrequency chosen at random from a probability 
distribution). Note that the interaction strength is symmetric between oscillators with 
i,j 7^ s, continuous in time and instantaneous. It depends on the natural or artificial systems 
that are modelled by Kuramoto oscillators, whether these features amount to appropriate 
idealizations or to an oversimplification. For example, the interactions between fireflies 
propagate with the very velocity of light. In view of their mutual distances on the trees 
the delay caused by the finite velocity of light should be completely negligible. When the 
interaction proceeds via pulses as in neural systems, the continuous interaction [as in system 
(H))] is also not adequate. 
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III. PHASE-LOCKED MOTION 



In this section we consider the conditions for having phase-locked motion, in which the 
phase differences between any pair of oscillators remain constant over time, after an initial 
relaxation time. In general two quantities are of interest for phase-locked motion, the fre- 
quency Q, in common to all oscillators in the phase-locked state, and the critical threshold 
in terms of the absolute value of the eigenfrequency of the pacemaker \uJs\ over the coupling 
strength K, i.e. \uJs/K\ij, that gives an upper bound for having phase-locked motion as a 
stable state. We analytically calculate (see the Appendix) Q for any dimension d of the lat- 
tice and for arbitrary range of the interaction (i.e. < a < oo). Imposing the phase-locked 
condition 

<^^ = Q , Vz = 0, ...,iV 
to system and using the fact that the sine is an odd function, one easily obtains 

where the r]iS are defined in Eq.Q. They contain the informations about the range of the 
interactions and about the topology on which the system is defined (i.e. rji is a function of 
a, it depends on the dimension d and on the boundary conditions of the lattice). Moreover, 
in case of periodic boundary conditions and in case of long-range interactions (a = 0), we 
find 

r], = 7] , V z = 0, . . . , iV , 

or, in other words, the normalization factor defined in Eq.(j21) becomes independent on i 
(i.e. the position on the lattice of the i-th oscillator). In this case we are able to prove that 
Eq.(jni) takes the explicit form 

The main result is that Eq.(0]) is independent on a. Moreover, in such cases the position 
of the pacemaker does not change the behavior of the system. Therefore we dropped the 
index s of the pacemaker in Eq.(^ . 

As a next step, we calculate the critical threshold \uJs/K\(j for a one-dimensional lattice 
and different range of interactions: analytically for next-neighbor {a — > oo) and for long- 
range {a = 0) interactions f section IlII AIIIII'C|) . numerically for intermediate coupling range 
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fsection [lIIE|) . Both quantities, the common frequency Q and the critical threshold \uJs/K\^, 
have an analytical generalization to higher dimensions d > 1 f section |IIIF|) . 
All numerical results (apart from those of Figure are obtained with an accuracy of 5 ■ 10~^ 
by integrating the system by the fourth-order Runge-Kutta method (with fixed time step 
At = 0.05). 

A. Kuramoto oscillators on a linear chain with next-neighbor interactions 

We consider (A^ + 1) Kuramoto oscillators coupled along a chain as indicated in Figure 
n with free boundary conditions at the positions and A^, and the pacemaker located at 
position s, < s < A^. For this simple topology, Eq.(j2)) takes the form 



where 6{x) is the Heaviside step function [6{x) = 1 for x > and 6{x) = for x < 0]. Eq.(jS)) 
inserted into Eq.Q gives the proper value of the common frequency Q as function of the 
position s and the eigenfrequency Ug of the pacemaker, the size of the system A^ as well as 
the range of the interaction, parameterized by a. 

To derive the critical frequency of the pacemaker in the case of next-neighbor interactions 
(a — ^ oo), we use the odd parity of the sine function. We derive a recursive relation 
between the phase differences of next-neighbor oscillators. Using the boundary conditions 
it is possible to express these recursive equations in terms of s, N, Q, and K. On the 
other hand, the sine function is bounded. These bounds induce the bound on the maximal 
pacemaker frequency \uJs/K\^ that is compatible with a phase-locked motion. During the 
recursions between phase differences we have to distinguish whether we move to the boundary 
to the right or to the left of the pacemaker s, placed between < s < A^; either the boundary 
conditions enter for oscillator A^ or oscillator 0. Thus if we move to the right, we obtain for 
the maximal pacemaker's frequency 
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2A^ - 2s - 1 



(6) 



so that all oscillators to the right (i?) of s approach a phase-locked state, and 



L 



K c 



s 



(iV-l)(l + e,) + l 
2s - 1 



(7) 



6 



as bound for oscillators on the left (L) of s to synchronize in a phase-locked motion. Since 
we are interested in a state with all oscillators of the chain being phase-entrained, we need 
the stronger condition given by 

(8) 
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For the pacemaker placed at the boundaries s = or s = we have 
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(9) 
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If the pacemaker is located in the middle of the chain, i.e. at s = N/2 for N even or 
s = (A^ ± l)/2 for odd, we have 



and 



K 



K 



c 



c 



(Ar-l)(l + 6jv/2) + l 

N -I 

(Ar-l)(l + e(iv±i)/2)- 
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if N is odd . 



(10) 



(11) 



Note that for = — 1 the critical ratio goes to zero for — > oo, while it approaches one 
for = and A^ — » oo, but the common frequency f2 goes to zero in this case for A^ — cxd. 
Moreover we see the stronger the coupling the larger \ujs\ may be up to which the system 
is able to follow the pacemaker. 



B. Ring topology and next-neighbor interactions 

For a oo and (A^+1) oscillators placed on a ring, it is easy to see that the normalization 
factor defined in Eq.Q is the same for all oscillators i [therefore Eq.(0]) gives the right value 
of n\. It can be written as 
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if A^ is odd 



(12) 



The derivation on the upper bound of the pacemaker's frequency \uj/K\^ for a ring 
proceeds in analogy to Eq.s (0), © and 0, apart from the fact that the boundary conditions 



are periodic 



The result is 
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Figure 2: Critical threshold Iw^/iTl^ as function of the position s of the pacemaker on a one- 
dimensional lattice (for both free and periodic boundary conditions). Here = 20, a = 100.0 and 
e<j = 0. Theoretical predictions [dotted, dashed and full lines from Eq.s®, and (^IJ) , respectively, 
and stars from Eq.@] are in perfect agreement with results obtain from numerical simulations 
(circles and squares). 

C. All-to-all coupling 

In case of all-to-all coupling [a = 0) the topology of the system has no longer an effect on 
the phase entrainment; moreover the normalization factor of definition ^ is rji = N , V i = 

0. . . . , A^, therefore the common frequency is given by Eq.(0]). Since we know that in the case 
of all-to-all couplings the mean- field approximation becomes exact, it is natural to introduce 
a quantity that usually serves as order parameter ^ for indicating the synchronized phase, 

1. e. 
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In our definition R e^^ differs from the usual order parameter only by the fact that the sum 
runs over j ^ s. In terms of the order parameter we can rewrite the system (Q) according to 

VL = u}s + {1 + es)KR sin {ijj - (ps) ^^^^ 
= ^ sin {(fs — ^i) + KR sin {i}) — Lpi) ,\fi^s 

As we see from Eq.s fll6p . all oscillators different from the pacemaker satisfy exactly the 
same equation. One possible solution is given by 

so that all phases of the oscillators different from the pacemaker are equal. Using system 
(|T6|l and the fact that the sine function is bounded, it is easily seen that the critical ratio 
\uj/K\^ is still given by Eq. ()14|) . valid for the ring and a oo. 



D. A switch to synchronization 

Let us summarize the results we obtained so far in Figure El The full lines represent the 
analytical results for \uJs/K\fj for the open chain with a oo, more precisely for a = 100, 
and the ring for a = 100, the circles and squares represent numerical data that reproduce 
the analytical predictions within the numerical accuracy. All results are obtained for = 20 
and eg = 0, they are plotted as a function of the pacemaker's position that only matters in 
case of the open chain and for a > 0. The horizontal line at \u!s/K\^ = {N + 1)/N obviously 
refers to the ring, the two branches (left and right), obtained for the chain, cross this line 
when the pacemaker is placed at s = N/2 in the middle of the chain. When the pacemaker 
is located at the boundaries s = and s = A^, we obtain two isolated data points, cf. the 
corresponding Eq.Q, close to the horizontal line. 

Let us imagine that for given A^ and the absolute value of the pacemaker's frequency \ujs\ 
and the coupling K are specified out of a range for which 
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such that the ratio is too large to allow for phase-locked motion on a chain, but small enough 
to allow the phase-entrainment on a ring. It is then the mere closure of the open chain 
to a ring that leads from non-synchronized to synchronized oscillators with phase-locked 
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motion. Therefore, for a whole range of ratios \uJs/K\^ no fine-tuning is needed to switch 
to a synchronized state, but just a simple change in topology. Because of its simplicity we 
suspect that this mechanism is realized in biological systems that need the possibility of an 
easy switch between synchronization and desynchronization. In our numerical integration 
we simulated such a switch and plot the phase portrait in Figure El The phase evolution (p as 
function of time is always projected to the interval [0, 27r): we use a thick black line for the 
phase of the pacemaker and thin dark-grey lines for the other oscillators. In the numerical 
simulation with T = 8000 integration steps over time we analyzed a one-dimensional lattice 
of Kuramoto oscillators with = 6, a = 5.0, s = 2, uOg/K = 1 and = 0. In the time 
interval from to "ON" {T/3) we see for the system motions with non-constant slopes: more 
steep for the pacemaker and the left part of the system (oscillators z = 0, 1) and less steep 
for the right part (oscillators i = 3,4, 5, 6). At the instant "ON" we close the chain, passing 
to a ring topology. The system almost instantaneously reaches a phase-locked motion (all 
phases moving with the same slope). At time "OFF" (2T/3) we open the ring, again, and 
the system shows qualitatively the same behavior as before the closure. 

Furthermore it should be noticed from Figure |21 that it is also favorable to put the pace- 
maker at the boundaries of an open chain to facilitate synchronization. For d = 1 the 
pacemaker then has to entrain only one rather than two nearest neighbors so that the range 
of allowed frequencies \u!s\ increases. 

^From these results it is natural to try to utilize the topology in a way that facilitates 
synchronization in artificial networks when synchronized states are needed. Given A^ oscil- 
lators and a number of po pacemakers, one can optimize the placement of the pacemakers 
when only a limited range of couplings is available and no other fine-tuning of parameters 
is feasible. 

E. Intermediate range of couplings 

Recall from section IIIII that the common frequency Q of phase-locked motion was inde- 
pendent on a in case of a ring and depending on a in a well-defined way [cf. Eq.sQ,©] 
in case of an open chain. The upper bound on the pacemaker's frequency, however, was 
analytically derived only for a; = and a ^ oo. For intermediate values of a we studied 
the behavior of \uJs/K\^ as a function of a numerically. The results are shown in Figures ID 
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Figure 3: Phase portrait of a system of Kuramoto oscillators on a one-dimensional lattice with 
N = 6, a = 5, s = 2, eg = and uJs/K = 1. The phase of the pacemaker is represented by 
the fat black line, while the phases of the other oscillators are drawn with thin dark-grey lines. 
We start from a chain (free boundary conditions), for which we are above the critical threshold 
given in Eq.®, here the stable state of the system is characterized by a non-constant frequency 
(the frequency of the oscillators are different from each other and vary during time). We observed 
the system for 8000 integration steps altogether. At the instant "ON" (T/3) we closed the chain, 
passing to a ring topology: now the ratio u)/K is less than the one given by Ea. H14() . so that the 
system can synchronize in the frequency domain (it reaches a phase-locked motion where all the 
phases have the same constant slope). At instant "OFF" {2T/2>) we opened the ring and the system 
falls again into a non-phase-locked motion, similarly to that seen in the interval [0,ON]. 

In case of the ring topology (Figure |3)) we observe a non-monotonic function of a with a 
minimum at um < 3 that seems to decrease as the system size increases, here from = 10 
to 30. The shape of \uJs/K\^ depends on the choice of the normalization with K/r]i in 
the interaction term of Eq.(P3). If we had chosen K/N instead of K/rji, we would find a 
monotonic decrease (inset of Figure |3)). Instead we see here that all-to-all and next-neighbor 
interactions lead to the same ratio of \uj/K\(j, while synchronization becomes more and more 
difficult for a around um- The existence of um is also supported by numerical results of 
for A^ = 100 oscillators, where it was found that the system approaches the same behavior 
as the original Kuramoto model with a = for a < 2 so that there is a different behavior 
for long- and intermediate-range couplings. We expect that our values for um decrease to 
ctM 2 for larger system sizes. 
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Figure 4: Critical ratio \{jjIK\q for a one-dimensional ring with e = as function of a for various 
system sizes, showing a non-monotonic behavior with a minimum at au ^ 3. The inset shows the 
same critical ratio for N = 10, 14, 20, but for the system with r] = N, chosen independently on 
a: now \uj/K\^ behaves as a monotonically decreasing function of a. 

Figure El shows the same critical ratio \uJs/K\(j as a function of the pacemaker's position s 
for an open chain, = 20, = 0, and various values of a. Here it is interesting to see that 
the shape of \u!s/K\(j drastically changes around the same value of aM < 3 as for the ring: 
for a < the critical fraction looks different iable for all positions s of the pacemaker 
along the chain, apart from the boundaries s = and s = N, while for a > it appears 
non-differentiable around s = N/2 and has a cusp. An analytical understanding of \ujs/K\q 
as a function of s and a is still missing. 
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s s s 

Figure 5: Critical ratio \uJs/K\(j as function of the pacemaker's position for an open chain and 
various values of a, indicating a qualitative change in the behaviour between a = 2.5 and a = 3. 

F. Kuramoto oscillators on regular structures in higher dimensions 

In order to generalize our previous results to regular structures in higher dimensions, 
we consider Cayley-trees and hypercubic lattices. Cayley-trees have z branches to nearest 
neighbors at each node, cf. Figure IHl apart from nodes in the outermost shell, where the 
number of nearest neighbors is only z — 1. Cayley-trees generalize linear chains from z = 2 
to z > 2 branchings at each node. They provide a geometry that often allows for exact 
analytical solutions, cf. the percolation problem on a Cayley-tree IHJJ. Cayley-trees are 
infinite-dimensional in the sense that the number of nodes in the boundary (outermost shell 
at distance R from the origin) are proportional to the volume, i.e. the total number of nodes 
in the tree. Throughout our calculations we place one pacemaker in the center of the tree 
and assign all other Kuramoto oscillators with eigenfrequency zero to the remaining nodes. 
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6 " m 

Pajek 

Figure 6: A Cayley tree with z = 3 branches and radius R = 5. 
We start from a modified version of the equations (QJ 

(Pi = SifiLJo^ "^Aji sin {(pj - (pi) ,Wi = 0,...,N , (17) 

where Aji is the adjacency matrix of a Cayley-tree. Naively we would expect that the 
larger the number of branches z and the larger the radius i?, the more difficult gets the 
synchronization of the tree. Therefore we calculate the common frequency and the critical 
threshold \jjQjK\^ as a function of z and R. It is convenient to assign the distance r from 
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the origin (center node) as coordinate of a node together with the path, along which it is 
reached from the origin, that is {r, si, S2, Sr) with si = l,..,z, Si = l,...,z — 1 for all 
shells i = 2, ...,r. Imposing the condition of phase-locked motion, we see that all oscillators 
at distance r satisfy the same equation, independently on the position within this shell 
(i.e. independently on the path along which they are reached from the center), so that 
we finally can drop si,...,Sr and keep only the distance, ipr,si,s2,...,sr = fr- Again one can 
derive recursive relations for phase differences between phases ip in neighboring shells (cf. 
the Appendix) and express the final difference between shell and shell 1 exclusively in 
terms of the parameters of the system, i.e. z and R, leading to 

n = uo + z J2{z - ly + {z - If-' j . (18) 

As it is seen from Eq. p8|) . Q goes to zero for i? ^ oo, 2; ^ oo or both. The critical 
threshold \ujq/K\^ is obtained from an expression for sin{(fr — fr-i) (cf. the Appendix) 
that is always negative and monotonically increasing as a function of r if c^o > (or always 
positive and monotonically decreasing as a function of r if c^o < 0) , so that it takes a 
minimum (maximum) at r = 1. It is then again the bound of the sine function that leads 
to the critical ratio 



K 



c 



\R-1 



(19) 



Similarly to the result for \ujs/ K\^ for the linear chain, this ratio goes to one for z — 5> 00 or 
-R — > 00 or both, while ^ in the same limiting cases. For the special case oi z = 2 the 
Cayley-tree reduces to a linear open chain with the pacemaker placed in the middle of the 
chain, and Eq. (fT^ reduces to Eq. lfTm) with = 0. 

Since we are interested in finite systems, we also consider Eq. fll9p as a condition on the 
critical size in terms of z and R which may not be exceeded for keeping phase-locked mo- 
tion if the pacemaker's frequency ujq and the coupling strength K are given. The fact that 
for sufficiently small z and R synchronization on a Cayley- tree is possible, is in qualita- 
tive agreement with a result, recently obtained for an ensemble of Rossler oscillators on 



a Cayley-tree 







Rossler oscillators that are individually even chaotic systems, also ap- 
proach a synchronized state on a Cayley-tree if z and R are sufficiently small. We can 

extend qualitatively all the results so far obtained for any dimension d, when system (P) is 
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Figure 7: Critical ratio 
a = 1.0, b) a = 100.0. 



uj^(s)y(s)/K for a two-dimensional lattice with = 6, 6^=0 and a) 
" c 



placed on a hypercubic lattice with {Nj + 1) oscillators in each dimension, so that we have 
an ensemble of THj^^^^Nj + 1) Kuramoto oscillators. The i-th oscillator's position is labelled 



by (i-dimensional vectors x^*^ with < a;^*'* < , V j = 1, . . . , ci. A phase-locked motion 



now occurs for 



J^)/K < {u (s) js)/K 



quency at position x^^^ = {x[^' , . . . ,xl^'). In Figure [3 we have sketched \uj^(s)y(s) / for a 
square lattice in two dimensions with the pacemaker at x^*-* 



c 



the critical ratio for the pacemaker's fre- 



(a;("), e [0,1,..., 6], 



y 



is) 



G [0, 1, ...,6], a = 1.0 (Figure[7ti) and a = 100.0 (Figure[7|D), = 0. The largest range 
for synchronized states is obtained if the pacemaker is either put in the center of the square 
(for open boundary conditions), or on a torus in two dimensions (as generalization of the 
closed ring in c? = 1 dimension). The different shape of the "hats" at the boundaries is due 
to the dependence of the critical ratio on a, in particular at the boundaries, as it was visible 
also in FigEl Recall that Eq. © was derived only for a — >^ oo. 



IV. ABOVE THE CRITICAL THRESHOLD 



In this section we study the phase evolution of the oscillators when the ratio of the 
pacemaker's eigenfrequency u over the coupling K lies above the threshold of phase-locked 
motion. As our numerical integration of the differential equations shows, there are still some 
regular remnants to phase-locked motion above the transition. Far above the transition the 
oscillators can no longer follow the pacemaker and get stuck to their eigenfrequency zero. 
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These features are manifest in a phase portrait ipi{t) of all oscillators i as a function of time, 
and in some average frequency as a function of the pacemaker's eigenfrequency uj; Q 
replaces the common frequency Q below the critical threshold. Its precise definition is given 
below. Here we consider only the one-dimensional ring-topology. We report our results in 
Figures |H1 and IHl 




t 




Figure 8: Phase portrait above the critical threshold {uj/K = 2 > \u/K\q) for = 6 oscillators 
on a ring with a = 5.0 and e = 0. The almost vertical black lines represent the pacemaker's phase, 
while the phases of the other oscillators are plotted with thin dark-grey lines, a) The pacemaker 
moves faster (steep slope) than the rest of the system that seems to move coherently with a common 
frequency; b) zooming into the phase portrait, it is clear that the phase of each of the oscillators 
different from the pacemaker oscillates around the common frequency with an amplitude inversely 
proportional to the distance from the pacemaker. 

Phase portrait 

Figure IHt^ shows the portrait of a non-synchronized ring with N = 6, a = 5.0, e = 0, 
and \uj/K\ = 2.0 > \uj/K\^. The phases ip are projected to the interval [0,27r). The fat 
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Figure 9: Dependence of the average frequency Q [Eq.(j221)] as function of the pacemaker's frequency 
uj. Numerical results refer to the case of a ring, with e = 0, liJ>0, i^ = l and to different sizes of 
the system (i.e. different values of N) and/or ranges of the interactions (i.e. different values of a). 
In this case, system is integrated using the eight-order Runge-Kutta method with variable time 
step. Each point is given by the average over twenty different trials (for which the initial phases are 
chosen at random from a uniform distribution) and the error bars represent the standard deviation 
from the mean values. Error bars appear only from = 50 on and for a > 0, because the value 
of to ~ 10'^ is not always large enough to ensure that the system reaches the stable state for the 
given initial condition. 



vertical lines represent the phase evolution of the pacemaker, while the thin dark-grey lines 
correspond to the phase evolution of all other oscillators which seem to almost coincide along 
these lines, but a zoom, shown in Figure IHb; displays their difference in amplitudes as well 
as the superposition of a fast and a slow motion. Rewriting the differential equation for the 
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i-th oscillator from system (Ql, with periodic boundary conditions, as 




sin {ipj - ipi) 



(20) 



we can distinguish the interaction between the i-th oscillator and the pacemaker from the 
interactions between the i-th oscillator and the rest of the system. Calculating from Eq. ()20|) 
the average frequency Vt of all the oscillators different from the pacemaker 



we see that Vt is given by a weighted average of the interactions between the pacemaker 
and the other oscillators, without any contribution from the mutual interactions between 
oscillators different from the pacemaker. 

Numerically we approximated the average frequency of Eq. ()2H) as 



with 10^ ~ ti ^ to ~ 10^. The values of ti and to are chosen sufficiently large to capture 
only the slow motion {2-n / (ti — to) ^ t^)- Note that reduces to when \uj/K\ < \uj/K\^, 
since in the phase-locked regime = , V i. 

Average frequency as a function of the pacemaker's eigenfrequency 

Moreover, it is of interest to study the absolute value of the average frequency as a 
function of the pacemaker's eigenfrequency uj, see Figure El where we have chosen K = 1 
for simplicity as well as u; > 0. For uj < ujc — 1 we obtain the linear increase of i7 = oc tu, 
as we expect from Eq.(^. For ^ oo we expect i7 ^ 0, since the oscillators can no 
longer follow the pacemaker when its eigenfrequency goes to infinity. However, we have no 
analytical understanding of the decay of Vl as an inverse power of tu, i7 oc cj^^ for uj ^ 1. 
The dotted envelopes mark f2 oc cu, while the dashed curves correspond to ^7 oc uj^^. The 
numerical results show a weak dependence on the interaction range (a varies between 
and 10), as well as a stronger dependence on the size of the system [for uj < uoc as 
predicted by Eq.(j3])], varying between 10 and 50. For the numerical integration of the 
differential equations we used the eight-order Runge-Kutta algorithm with adapted time 




(21) 




(22) 
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steps in order to accelerate the integration with variable time step. This introduces some 
"noise" in the numerical estimate of the first derivatives. Each data point corresponds to an 
unweighted average over twenty different choices of the initial phases (i.e. chosen at random 
from a uniform distribution). (To achieve a complete independence on the choice of initial 
conditions and a full stabilization would need a longer simulation time.) The size of the 
errors increases with as indicated in Figure IHl 

V. SUMMARY AND CONCLUSIONS 

We considered ensembles of Kuramoto oscillators with zero eigenfrequency apart from 
one, the so-called pacemaker, with eigenfrequency larger than zero and coupled to the oth- 
ers in a symmetric (e^ = 0) or asymmetric (— 1 < e < 0) way. The oscillators were assigned 
to d-dimensional regular topologies, in particular to a one-dimensional open chain, a one- 
dimensional ring, a Cayley-tree with z branches and a hypercubic lattice in d dimensions. In 
general, we studied the conditions for having phase-locked motion as a function of the system 
size iV, the coupling strength the pacemaker's eigenfrequency lo^ and the asymmetry pa- 
rameter eg. The interaction range between the oscillators was varied between next-neighbor 
and all-to-all couplings via the parameter a. We derived the common frequency of the 
phase-locked motion, and the critical threshold \iOslK\^^ analytically for the limiting cases 
a oo and a = 0, in particular for a chain with arbitrary position of the pacemaker and 
a ring in one dimension, as well as for a Cayley-tree with coordination number ^ > 2. In- 
termediate coupling ranges < a < oo and the phase evolution above the critical threshold 
for phase-locked motion were treated numerically. 

In the large limit, phase-locked motion is impossible in all cases we have considered, 
that is, for both chain- and ring- topologies, and for symmetric (e^ = 0) and asymmetric 
(— 1 < ts < 0) couplings of the pacemaker. For (e = —1), the common frequency Q con- 
verges to the pacemaker's Ug, but \uJs/K\^ goes to zero for ^ oo, whereas ioi e = fl 
goes to zero while \uJs/K\^ goes to one for oo. Similarly, on a Cayley-tree, which we 

studied for = 0, Q and \ujq/K\^ go to zero for N ^ oo, independently of whether the 
limit is realized as z — oo, or i? — > oo, or both. However, we were not only interested in 
the "thermodynamic" limit and do not mean the transition to the desynchronized phase 
for finite A^ in the thermodynamic sense of a phase transition. Our main result for finite 
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and — 1 < e < was that for a whole range of ratios \uJs/K\ it is possible to induce 
synchronization by a mere closure of a chain to a ring. (For ratios below this range the 
system synchronizes even for an open chain with pacemaker at position s, for ratios above 
this range synchronization is neither possible on a ring.) 

Two interesting results were obtained from the numerical integration of the differential 
equations. Firstly, as a function of the interaction range (parameterized by a), the critical 
ratio \uj / K\(j showed a minimum at some intermediate value around ctAf < 3 (for given 

and e that we considered). This non-monotonic shape was obtained for our choice of 
normalization r^j, defined in Eq.Q. For the pacemaker it is therefore as easy to entrain the 
phases of the other oscillators if it only interacts with nearest neighbors as if it interacts 
with all other oscillators as nearest neighbors (mean- field case) . This is well understandable 
from Eq.s ()24|) and (^Sjl in the Appendix, from which it is seen that it is sufficient for the 
pacemaker to entrain its nearest neighbors in order to entrain all other oscillators. Formulas 
fl24|l and (j25|) were derived under the condition of nearest-neighbor interactions. Therefore 
the argument does not go through for intermediate values of a, for which the numerical 
results show that the pacemaker has a "harder job" to entrain the remaining ensemble. 

Secondly, above the transition from the phase-locked motion to the desynchronized phase, 
we found some regular remnants of phase-locked motion. It is still the pacemaker that 
determines the fast frequency of phase fluctuations of the others around some average value Vt 
(also determined by the pacemaker), and it is the distance to the pacemaker that determines 
the amplitude of these fluctuations. The value of the average frequency f2 varies slowly as 
compared to the pacemaker's frequency uj. As expected, VL approaches zero for \uj/K\(j oo 
as an expression of the fact that the system can no longer follow the pacemaker in this limit. 
The simple power-law of the decrease of VL remains to be understood analytically. 
The sensitive dependence of synchronization on the topology in a certain range of parameters 
may be exploited in artificial networks and is -very likely- already utilized in natural systems, 
in which a switch to a synchronized state should be easily feasible. 
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APPENDIX 

Kuramoto Model on a d-dimensional lattice 

We start with the system given by Eq.s(IT)) and defined on a dimensional lattice. 
Imposing the phase-locked condition 



we have 



where 



ip^ = n ,Wi = 0,...,N 



K 

Q = Si^s^^s + (1 + Si,s^s) — , V 2 = 0, . . . , 

Vi 



sin {ipj - 



Common frequency in the phase-locked regime 

As it is easily seen from the odd parity of the sine function, 

f2j = and Vti = , 

SO that 
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while for j 7^ s 

=^ # E,v. E.^, r/. = ^ + (1 + e.iV)^]. 

which imphes 

The common frequency then follows as in Eq.(jn}. 

When the topology has a symmetry such that the normalization factor defined in Eq.(j21) is 
independent of the index i 

r], = 7] , V z = 0, . . . , iV , 

then 

(1 + eN)^ri + ^^7] = {I + eN) Nt] + N^r] = Nr] {N + eN + 1) , 
thus we obtain Eq.Q. 

Such symmetries are realized for a = 0, and for arbitrary a if we have periodic boundary 
conditions. 



Critical ratio \lOs/K\(-, for the one-dimensional lattice 



In the case of next-neighbor interactions (i.e. for a — > 00), on an open chain, we can 
solve the system (P) recursively. We start from the equation for the pacemaker (assuming 
< s < and 7^ —1) and move first to the right of the pacemaker 



n = ^ 

" 2 



sin {ifis - ^s+i) + sin {ips+2 - Vs+i) 
4^ sin {lPs+2 - ^s+i) = ^ + sin (v^^+i - cp, 
For all 1 < J < — s it then follows that 



n = f 



sin (<^s+(j-2) - V5s+(j_i)) + 
+ sin {cps+j - 
^ sin {(ps+j - ips+(j-i)) = 
= ^ + sin (y^s+ij-i) - Vs+ij-2)) = 
= U - l)x + sin {(Ps+i~^s) 
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In particular for j = N — s, we have 

2Q 

sin {ipN - (pN^i) = (A^ - s - 1) — sin {ips+i - ips 

K 

Imposing the boundary conditions 

Q = K sin {ipN-i — ^n) 

K ' 



<^ sin {ipN-i - V^iv) - ^ 



we obtain for all 1 < j < — s 

sin(<^,+,-<^,+(,_i)) = (2s-2iV-l)^ + 2j|: , (23) 

which implies Eq.(jn)) by imposing the condition that the sine is a bounded function (| sinx| < 
1). When ujs > Q [uJs < 0], Eq. ()23p is always negative [positive] and has its minimum 
[maximum] value for j = 1. So if 



\sm{ips+i-ips)\<l , V 07,^0 (24) 



sin {ips+i - ^s) > -1 , if > 
sin (vj^+i - v^s) < 1 , if < 
then 

sm(</^«+j -^s+{i-i)) > -1 ,if^«>0 

I sm [ips+j - (fs+{j-i)) I < 1 , V 7^ , 

sin {ips+-j - V^s+o-i)) < 1 , if o;^ < 

(25) 

Vl<j<A^ — s. In other words, if the first neighbor on the right hand side of the 
pacemaker (the (s + l)-th oscillator) follows the pacemaker [i.e. if Eq. ()24|) is satisfied], then 
all the others follow the pacemaker with the same velocity [i.e. Eq. ()25|) is also satisfied]. 
In the same way we obtain Eq.© for oscillators to the left of the pacemaker. Of course we 
need that both the left and the right parts are synchronized in frequency, hence Eq.® is 
the global solution. 

Eq.®, for which the pacemaker is at the first (s = 0) [or at the last (s = A^)] position of 
the chain, is obtained in the same way, but in these particular cases only a motion to the 
right [or to the left] is possible because the left [or the right] part does not exist. In case of 
the ring the procedure for deriving Eq.(^H) is the same, apart from applying the different 
boundary conditions, periodic ones in this case. 

For a = and phase-locked conditions, the system ((T]) and Q becomes 

= Si^sUJs + (1 + (^i.sCs) ^ sin (v?j - (fi) . 
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Introducing the definition of Eq. ()15|) (whicli serves as an order parameter in a mean- field 
approximation that here becomes exact), we write 

Q = Us + {I + es)KRsm {ip — ips) for i = s 
i7 = ^ sin {(ps - (fi) + KR sin {ip - ipi) , \/ i ^ s 

Here all oscillators different from the pacemaker satisfy exactly the same equation. One 
possible solution is that all the phases of the oscillators different from the pacemaker are 
equal 

ipi=ij , W i s , 

from which it follows that 

K 

f2 = — sin -Vi) , V 2 7^ s . 

Using the last equation, the fact that the sine-function is bounded and finally the value 
of the common frequency, previously calculated in Eq.(j3|), it is straightforward to derive 
Eq.dHD. 



Kuramoto Model on a Cayley Tree 



Consider a node at distance R from the pacemaker, reached along the path 



(si, S2, S3, ... , sr_i, sr) from the central node. Here Si = 1, . . . , z , while Sj = 1, 
\/i = 2, . . . , R. For this node we can write Eq. (fTTj) as 

Q = 



,z -1 , 



K 



= Ksin {(pR-l,si,S2,...,SR.i - 'fR,si,S2,...,SR.uSR) 
^ sin{Lpji_i^s^^S2,...,sii_i - ^R,si,S2,...,sn_i,SR) 

from which we can see that all oscillators at distance R, independently on the path along 
which they are reached from the center, satisfy the same equation. 
Proceeding in an analogous way to before, we can write 

Q = 

,S2,...,SR-1,SR V^i?-l,Sl,S2,---,SJS-l) + 

+ sin {iPR-2,Sl,S2,...,SR-2 ~ fR-l,Sl,S2,...,SR-l, 



'fR-l,Sl,S2,---,SR-2,SR-l 



K_ 

z 



-^)^ + sin {^R-2,si,S2,...,SR^2 ~ '■PR-1,suS2,...,sr^j, 



^ sin (v9jj_2,si,S2,...,SH_2 - V^i?-l,si,S2,...,sj 



R.,) = ^[Z + {Z - 1)] , 
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next 



V5i?-2,si,S2,...,s_R_ 



K 



Z]L-i=l sin(v?/?-l,si,S2,...,Sfl_2,Sfl_i - V^R-2,si,S2,...,Sfl-2) + 
+ sin(v9/j_3,si,S2v,«ii-3 ~ V'fi-2,si,S2,...,Sfl_2^ 



K 



{Z-1)^[Z+{Z- 1)] + sin (</?K-3,Si,S2,...,«fl-3 - </'iJ-2,.l,.2,...,Sfl-3,^fl-2) 
sin (<^i?-3,.i,«2,...,.K-3 - <^iJ-2,.l,.2,...,.«-3,^fi-2) = + ^(^ - 1) + - 1)^] 



and so on, until we arrive at 



sin {ipa^ 



(t+l),si,S2v.SH_(t+i) V'fi-t,si,S2,...,Sfl_(t4-i),Sfl_ty 



t-i 



g=0 



For t = i? — 1, we have 



R-2 



sin((^o - y^Mi) = i ^ 



but also, at the center of the Cayley tree, 



g=0 



iJ-2 



f] = 00 = ^0 + ^ sin((/.i,,, -v^o)=^o-^^ + (■2-1)''"'] • (26) 

si=l 9=0 

Eq. (fTHjl is easily obtained from Eq. (P^ . 

As we have seen, all oscillators at the same distance r from the pacemaker satisfy the same 
equation. We can write 



sin {ipr - ipr-l) 



K I + z E,="o (^-1)^ + (^-1)^-^ 
suppressing the path that was followed to reach the node. It is convenient to rewrite this 
equation according to 

\R~r 



Z 22g=0 


{z-iy + {z-iy 


1 






1 




- (z-iy-i 


1 






1 


1 + zJ2f=o (^- 





\R-1 



js=0 



\R-1 



\R-1 



^r-2 



for R>2. Finally we obtain 

UJo 



sin {(pr - (fr-l] 



1 + z e:=o (^ - 1)^ 



K{z- 1)^-1 L 1 + ^ ^fj _ 1). + _ i)R^i 



- 1 



(27) 
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When ujo > 0, Eq. ()27p is always negative, monotone and increasing and takes its minimum 
value for r = 1, while when < 0, Eq. (jTfj) is always positive, monotone and decreasing 
and takes its maximum value for r = 1, namely 

1 



sin {(pi - v^o) 



K 



In order to derive the maximal value for \ujq/K\, we use the bound on sine function, from 
which we find the critical fraction of Eq. ()19|) . 
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